Relaxation of finite perturbations: beyond the Fluctuation-Response relation 
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We study the response of dynamical systems to finite am- 
plitude perturbation. A generalized Fluctuation-Response re- 
lation is derived, which links the average relaxation toward 
equilibrium to the invariant measure of the system and points 
out the relevance of the amplitude of the initial perturbation. 
Numerical computations on systems with many characteris- 
tic times show the relevance of the above relation in realistic 
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LEAD PARAGRAPH 

Understanding the behavior of a dynamical sys- 
tem out of its equilibrium is a crucial issue of sta- 
tistical physics. In the case of an infinitesimal 
perturbation that shifts the system out of equi- 
librium, the classical Fluctuation-Response the- 
orem allows to determine the linear response of 
the system in term of its equilibrium properties, 
i.e. correlation functions. While the behavior 
of infinitesimal perturbations gives relevant infor- 
mation for problems of statistical mechanics, for 
climate and geophysical models the main goal is 
to characterize the relaxation of large perturba- 
tions, which can not be obtained from the linear 
response theorem. We present here a general- 
ization of Fluctuation-Response relation, which 
holds for finite amplitude perturbations, provid- 
ing a tool for extracting non equilibrium behav- 
ior out of equilibrium features of the system. We 
also discuss the non trivial role of the amplitude 
of perturbations in systems where many charac- 
teristic time scales are present. 



I. INTRODUCTION 

The Fluctuation- Response (F/R) relation has a deep 
relevance in statistical physics and more generally in sys- 
tems with chaotic dynamics (in particular in hydrody- 
namics [1] ) . The relevance of a connection between "non 
equilibrium" features (i.e. response to an external per- 
turbation) and "equilibrium" properties (i.e. time cor- 
relations computed according to the invariant measure) 



is well known in statistical mechanics. We can mention 
the important Green - Kubo formulas in the linear re- 
sponse theory [2]. Beyond statistical physics, an other 
field where the F/R problem has an obvious relevance 
is climate research [3]. One of the key problems is the 
possibility to understand the response of the present cli- 
mate to some violent changes (e.g. a volcanic eruption). 
The essential point is the possibility that the recovery of 
the climate system from a perturbation (response) can 
be estimated from its time history (correlations time of 
the unperturbed system). 

Assuming that the system is mixing and has invari- 
ant probability density function (pdf) p(x), it is possible 
to derive the following F/R relation. Let us denote by 
x(t) = (xi(t), ■ ■ ■ , Xiv(t)) the state of the system at time 
t. If at the initial time t = the system is perturbed by 
Sx(0) = (Sxi(0), ■ ■ ■ , &£jv(0)), the average evolution of 
the perturbation (6xi(t)} with respect the unperturbed 
trajectory is 



where 



, Sxj(t) , 



( Xi (t)f 3 (x(0))} 



and the function fj depends on p(x) as 



dlnp(x) 

dxj 



(1) 



(2) 



(3) 



In Section II we will give a complete derivation of the 
above formulas. 

As far as we know, the F/R problem had been studied 
only for infinitesimal perturbations. For statistical me- 
chanics problems it is relevant to deal with infinitesimal 
perturbations on the microscopic variables. In a simi- 
lar way this problem has importance in many analytical 
approaches to the statistical description of hydrodynam- 
ics where Green functions are naturally involved both in 
perturbative theory and closure schemes [1,4]. 

On the other hand in geophysical or climate problems 
the interest for infinitesimal perturbation seems to be 
rather academic, while the interesting problem is the be- 
havior of relaxation of large fluctuations in the system 
due to fast changes of the parameters. 
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In this paper we want to address the problem of the 
F/R relation for non infinitesimal perturbations. In Sec- 
tion II we will show that it is possible to generalize the 
F/R relation to large perturbations, involving rare events 
of the invariant measure. Section III is devoted to a dis- 
cussion on the connections, and differences, between our 
approach and well known results in dynamical system 
theory. In Section IV we will discuss the application 
to systems involving a single characteristic time, while 
Section V is devoted to system with many characteris- 
tic times. Section VI is devoted to conclusions and the 
Appendix VII contains some technical remarks. 



where 



F(x Q ,5x ) 



p(x Q - Sx ) - p(x ) 



P0o) 



(8) 



For an infinitesimal per- 

turbation 6x(0) = (Sxi(0) ■ ■ ■ fejv(O)) expanding (8) to 
first order one ends with the expression 



(5xi(t)) =-$>i(t) 



t=o 



(9) 



II. THEORETICAL BACKGROUND 

In the following we will consider a dynamical system 
with evolution x(t) = <p t x(Q) of the iV-dimensional vec- 
tor x. For generality, we will explicitly consider the case 
in which the time evolution can also be not completely 
deterministic (e.g. stochastic differential equations). We 
will assume the existence of an invariant probability dis- 
tribution p(x) and the ergodicity of the system so that 



lim 



A(x(t))dt = p{A) 



A(x)p(x)dx 
(4) 



for any (smooth enough) observable A. 

Our aim is the understanding of the mean response 
(8A(t)) of a generic observable A initially perturbed with 
5A{0). The first step is the study of one component of 
x i.e. (Sxi(t)) with an initial (non-random) perturbation 
Sx(0) = 5xo. Introducing the probability of transition 
from (jeo, 0) to (x, t), W(x , — > x, t) (for a deterministic 
system we have W(xo,0 — > x,t) = 5{x(t) — 0*xo)), we 
can easily write an expression for the mean value of the 
variable computed along the perturbed trajectory x'^t) — 
Xi(t) + Sxi(t): 



J J Xip (xo)W(xo, — > x, t) dxdx 



(5) 



where p'(x) is the initial distribution of perturbed sys- 
tem, which is related to the invariant distribution by 
p'(xo) = p(xo — Sxo)- Noting that the mean value of 
Xi{t) can be written in a similar way: 

(xi(t)} = j j Xip(x o )W(xo,0 -> x, t) dx dx a (6) 



one has: 



(Sxi(t)) = J J x l 



p(x - Sx ) - p(x Q ) 



P(xo) 

x W(xo, — ► x, t) dx dxo 
= (xi(t)F(x ,5x )} 



(7) 



which defines the linear response 

dlnp(x) 



Rij(t) = -(xi(t) 



dx. 



) 



(10) 



t=o 



of the variable x% with respect to a perturbation of Xj. 
Relation (10) is the generalization for non Hamiltonian 
systems of the well known fluctuation/response (F/R) 
relation [2]. 

Let us note that in the general case the invariant mea- 
sure p(x) is not known, so the equation (10) gives just 
a qualitative information. In the case of Gaussian distri- 
bution, p(x) factorizes and the linear response recovers 
the correlator 



(x l (t)x 3 (0)) - {xj){ Xj ) 
(xjXj) — (xj) (xj) 



(11) 



In the case of finite perturbations, the F/R relation (7) 
is typically non-linear in the perturbation 5xq and thus 
no simple relations analogous to (10) exist. Nevertheless 
we can disentangle the different contributions in the re- 
sponse (7) by studying an initial perturbation whose only 
non-zero component is the j — th one, 



5 {]) x{Q) = (0,---,0,<5x 3 (0),0, •••,()) 



(12) 



We therefore generalize the F/R relation (10) to non- 
linear response of Xi to a perturbation on the j variable 
as 



RiAt) = (Xi{t)f j{0)) 



where fj is given by 



p(x o -d^x(0))-p(x o ) 
fi{ 0) ~ p(x )S Xj (0) 



(13) 



(14) 



The explicit prediction of the response from (13) re- 
quires the analytic expression of the invariant pdf, which 
is in general not known. Nevertheless (7) guarantees the 
existence of a link between equilibrium properties of the 
system and the response to finite perturbations. This 
fact has a relevant consequence for systems with one sin- 
gle characteristic time: a generic correlation (e.g. the 
correlation (11)) in principle gives informations on the 
relaxation time of finite size perturbations, even when 
the invariant measure p is not known [5]. 
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III. REMARKS ON THE CONNECTIONS 
BETWEEN F/R RELATION, DYNAMICAL 
SYSTEM THEORY AND STATISTICAL 
MECHANICS 

Since the F/R relation involves the evolution of dif- 
ferences between variables computed on two different re- 
alizations of the system, it is natural to conclude that 
this issue is closely related to the predictability problem 
and, more in general, to chaotic behavior. Actually, a 
detailed analysis shows that the two problems, i.e. F/R 
relation and predictability, have only a very weak con- 
nection. For the sake of completeness, we shortly discuss 
here the analogies and differences between these two is- 
sues. 

The typical problem for the characterization of pre- 
dictability is the evolution of the trajectory difference 
Sx(t), in particular of (In |&e(£)|) which defines the lead- 
ing Lyapunov exponent A. For small |<5a?(0)| and large 
enough t one has 



<ln|<5as(t)|> ~ln|<fcc(0)|Ai 



(15) 



On the other hand, in F/R issue one deals with averages 
of quantities with sign, such as (Sx(t)). This apparently 
marginal difference is very important and it is at the ba- 
sis of the famous objection by van Kampen related to the 
standard derivation of the linear response theory [6]. In 
a nutshell, using the modern dynamical systems termi- 
nology the van Kampen's argument is as follows. Since 
in presence of chaos |5aj(t)| grows exponentially in time, 
it is not possible to linearize (8) for time larger than 
(1/A) ln(A/|&c(0)|), where A is the typical fluctuation 
of the variable x. As a consequence, the linear response 
theory is expected to be valid only for extremely small 
and unphysical perturbations, in clear disagreement with 
the experience. A solution of this apparent paradox was 
proposed by Kubo who suggested that "instability [of 
the trajectories] instead favors the stability of distribu- 
tion functions, working as the cause of the mixing" [7]. 
More recent works have demonstrated the constructive 
role of chaos in F/R relation and the non relevance of van 
Kampen's criticism [8,9]. The objection by van Kampen 
remains nevertheless relevant for numerical computations 
of F/R relation (see Appendix). 

Fluctuation/response relation was developed in the 
context of statistical mechanics of Hamiltonian systems, 
but it also holds for non conservative systems, and even 
non deterministic systems (e.g. Langevin equations) and 
has no general relation with "chaotic quantities" such as 
Lyapunov exponents or Kolmogorov-Sinai entropy This 
generated in the past some confusion about the appli- 
cability of F/R/ relation. For example, some authors 
claimed (with qualitative arguments) that in fully devel- 
oped turbulence there is no relation between equilibrium 
fluctuations and relaxation to equilibrium [10] while the 
correct statement concerns the non validity of the sim- 
plified relation (12) which holds only for systems with 
Gaussian statistics. 



Thanks to its general validity and robustness, the F/R 
relation has been also used to obtain informations on 
the unknown invariant measure p{x) on the basis of the 
linear response Rij(t). An important example comes 
from the field of disordered systems where the F/R had 
been applied to the study of aging phenomena [11]. 

Concluding this short discussion on the connections be- 
tween F/R relation, dynamical system theory and statis- 
tical mechanics, we mention recent results about rigorous 
derivation of the Onsager reciprocity relations [12] and 
the macroscopic fluctuation theory for stationary non- 
equilibrium states [13] in a class of stochastic models de- 
scribing interacting particles systems. 



IV. SYSTEMS WITH A SINGLE 
CHARACTERISTIC TIME 

Let us start by studying two examples of systems with 
a single characteristic time: a deterministic chaotic sys- 
tem (the Lorenz model) and a nonlinear Langevin pro- 
cess. 

We first consider the Lorenz model [14] 

dx , 

dy 



— = —xz + rx — y 
dt y 

dz 

— — xy — bz 
dt y 



(16) 



with standard parameters for chaotic behavior: b = 8/3, 
<t = 10 and r = 28. The correlation function (11) for the 
variable z, shown in Fig. 1, qualitatively reproduces the 
behavior of the response to different sizes of the pertur- 
bation of the z variable, ranging from infinitesimal ones 
up to the size of the attractor. The accuracy does not 
increase when decreasing the perturbation because the 
invariant distribution is not Gaussian (see Fig. 1) and 
thus the general correlation (10) should be used. We ob- 
serve that the use of (10) instead of (11) is in general 
much more difficult because the invariant distribution is 
in general non factorable. 

To better illustrate this point, let us now consider a 
system whose invariant probability distribution is known. 
In this case we can quantitatively compare the differences 
between the responses to infinitesimal and finite pertur- 
bations. Our example is provided by the stochastic pro- 
cess x(t) determined by 



dx 
~dt 



dU(x) 



(17) 



where £(t) is a white noise, i.e. a Gaussian process with 
(£(*)) = and (£(*)£(*')> = S(t - t'). The invariant 
probability distribution is [15]: 



p(x) = Ne- U ^/ D 



(18) 
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where Af is fixed by normalization. 

A Gaussian pdf is obtained using U(x) = x 2 /2 which 
corresponds to t he linear Ornstein-Uhlenbeck process 
dx/dt = — x + >/2D^(t). Our example uses a modified 
version of the Gaussian case, 

U= { ff x '\ X \ < \ (19) 
\x\ — 5 , |a;| > 1. 

The resulting pdf, shown in the inset of Fig. 2, has 
a Gaussian core, with exponential tails. Figure 2 also 
shows the response function for an infinitesimal and for 
a finite size perturbation. For both perturbations, the 
response function measured from the perturbed trajecto- 
ries is exactly predicted by statistics of the unperturbed 
system according to (13), while the Gaussian correlation 
C{t) = (x(t)x(0))/a 2 gives only an estimate of the relax- 
ation time. By construction, the pdf of this system has 
larger tails than in the Gaussian case, thus large fluctu- 
ations decay slower than small ones. In the linear case 
the mean response is simply R(t) — exp(— t) and does 
not depend on the amplitude of the initial perturbation 
6x(0). 

The results obtained for the Lorenz model and for the 
nonlinear Langcvin equations suggest that if only one 
characteristic time is present, the existence of the F/R 
relation allows for some qualitative results even in the 
absence of precise knowledge of p, both for infinitesimal 
and finite perturbation. 

V. SYSTEMS WITH MANY CHARACTERISTIC 
TIMES 

In systems with many characteristic times, different 
correlation functions do not show the same behavior, i.e. 
depending on the observable one can observe very dif- 
ferent time scales, corresponding to the different decay 
times of the correlation functions Cjj — (xj(t)xj(0)} [5]. 
In addition, at variance with systems with one single time 
scale, here the amplitude of the perturbation can play a 
major role in determining the response, because different 
amplitudes may affect features with different time prop- 
erties. 

The link between equilibrium and relaxation proper- 
ties established by the F/R relation (13) suggests that 
it is possible to relate different relaxation rates with the 
time scales measured by means of correlations. Consider 
the case of an observable A which depends on all the 
variables of the system {x\, ■ ■ ■ , Xn}- For infinitesimal 
perturbations, a straightforward generalization of (1,2) 
gives: 

(SA(t))=J2(A(x(t))f j (x(Q)))Sx j (0) (20) 

In the case of finite perturbations, as stressed in Sect. 
(II), it is possible to write a F/R relation: 

(SA(t))=(A(x(t))F(x(0),Sx(0))) (21) 



in which, at variance with (20), all the variables are 
mixed. In (21) the relaxation properties depend explic- 
itly on the initial perturbation Sx (0). 

Depending on the choice of A(x), different perturba- 
tions on A correspond to different amplitudes of the per- 
turbations on each variable Xj. Consequently, one can 
think that it is possible to associate each perturbation to 
a certain subset of variables which are mainly perturbed. 
The relaxation of (SA(t)) will be ruled by the character- 
istic time of that particular subset. 

In order to illustrate this issue we consider a shell 
model for turbulence [16]. Shell models are a simplified 
model for turbulent energy cascade, that describe the dy- 
namics of velocity fluctuations at a certain scale l n = fc" 1 
with a single shell- variable u n . Wave-numbers k n are geo- 
metrically spaced as k n = fc A™, allowing to cover a large 
range of scales with relatively few variables. A quadratic 
interaction between neighbor shell reproduces the main 
features of three-dimensional turbulence. The specific 
model we will use is 

{jtt + ^™) Un = i [ kn+lU n+i Un + 2 ~ e k n u n+1 ul L _ 1 

+ (1 - e)k n -iU n -2U n -i] + f n (22) 

where v is the molecular viscosity, /„ is an external forc- 
ing which injects energy at large scale, and e is a free pa- 
rameter. In order to have the correct conservation laws 
(energy and helicity) in the inviscid unforced case one 
has to fix e = 1/2. The observable considered is the to- 
tal energy E(t) — \ J2n=i \ u n(t)\ 2 which is the conserved 
quantity in the inviscid, unforced limit [16]. 

In order to study the response to perturbations with 
different amplitude on E, we consider the following per- 
turbed systems labeled with i = 1,..,N: u„\t) = 
u n (t) + Sun\t) where the initial perturbations Sun\o) 
are set in the following way: 

r (i), n > f , 1 < n < i - 1 

^ (0) = \v/(KF> ,*<n<N W 

This corresponds to a set of initial perturbations of the 
energy 

1 N 

<^(0)) = -]T(K| 2 } (24) 

n—i 

Such a perturbation is motivated by the fact that in 
the unperturbed system the energy is distributed among 
the shells according to the Kolmogorov scaling (|wn| 2 ) ~ 
kn 2 ^ , and the smaller scales give smaller contributions 
to the energy E(t). Thus it is natural to assume that a 
small perturbation of the energy will affect mainly the 
small scales. 

For each perturbation SEi, the average response of en- 
ergy 
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W)/ \Et 1 l^ ) (0)P-K(0)l 2 / 

reveals a close relation with the time correlation of the 
corresponding largest perturbed shell Ui(t), as shown in 
Fig. 3. A measure of the relaxation time can be provided 
by the halving times Ti/ 2 of the mean response, at which 
(5Ei(T 1/2 )) = l/2{5Ei(0)}. The dependence of response 
times on the amplitude of the initial perturbation, shown 
in Fig. 4, reflects Kolmogorov scaling for characteristic 
times t„ <~ k n 2 ^ 3 <~ u 2 n <~ 5E n 

Ti /2 - . (26) 

The above results on the shell model show that the 
response to a finite size perturbation of a system with 
many characteristic times may depend on the amplitude 
of the perturbation. Thanks to the existence of F/R 
relation is possible to establish a link between relaxation 
times of different perturbation and characteristic times 
of the system. 

VI. CONCLUSIONS 

Starting from the seminal works of Leith [3,17], who 
proposed the use of F/R relation for understanding the 
response of the climatic system to changes in the exter- 
nal forcing, many authors tried to apply this relation to 
different geophysical problems, ranging from simplified 
models [18], to general circulation models [19,20] and to 
the covariance of satellite radiance spectra [21]. In most 
of the applications it has not been taken into account 
the limits of applicability of the F/R relation which has 
been used as a kind of approximation. We have shown 
that a F/R relation holds under very general conditions. 
The derivation in Section II clearly shows the limits of 
applicability in its simplest form (i.e. the Gaussian ap- 
proximation (11)). 

Our main result is the demonstration that an exact 
fluctuation/response relation holds also for non infinites- 
imal perturbation. This relation involves the detailed 
form of the invariant probability distribution. In par- 
ticular, in order to predict the mean response to large 
perturbations, one needs a precise knowledge of the tails 
of the pdf. 

We believe that this generalization of the usual linear 
response theory can be relevant in many applications. As 
an example, we can mention climate research, where our 
results imply the possibility, at least in principle, to un- 
derstand the behavior of the system after a large impul- 
sive perturbation (e.g. a volcanic eruption) in terms of 
the knowledge obtained from its time history. Of course 
one has to take into account the strong limitations due 
to the need to have a good statistics of rare events. 



ACKNOWLEDGMENTS 

This work has partially supported by MIUR (Cofi- 
nanziamento Fisica Statistica di Sistemi Complessi Clas- 
sici e Quantistici) . We acknowledge the allocation of 
computer resources from INFM Progetto Calcolo Par- 
allelo. Angclo Vulpiani acknowledges support from the 
INFM Center for Statistical Mechanics and Complexity 
(SMC). We thank M. Falcioni for useful remarks. 

VII. APPENDIX 

In this appendix we want to discuss how van Kam- 
pen criticism is relevant for the numerical evaluation of 
infinitesimal response function. 

In numerical simulations, Rij(t) is computed perturb- 
ing the variable x% at time t — to with a small pertur- 
bation of amplitude 5xt (0) and then evaluating the sepa- 
ration 8xi{t) between the two trajectories x(t) and x'(t) 
which arc integrated up to a prescribed time t\ = to + At. 
At time t — t\ the variable Xi of the reference trajectory 
is again perturbed with the same Sxi(0), and a new sam- 
ple Sx(t) is computed and so forth. The procedure is 
repeated M 3> 1 times and the mean response is then 
evaluated according to (11). 

In presence of chaos, the two trajectories x(t) and 
x'(t) typically separate exponentially in time and the 
perturbed system relaxes to the unperturbed one only 
in average, therefore the mean response is the result of 
a delicate balance of terms which grow in time in differ- 
ent directions. The average error in the computation of 
Ri.jit) typically increases in time as e i( ^*/ 2 j\[M \ where 
L(2) is the generalized Lyapunov exponent [16]. Thus 
very high statistics is needed in order to compute B+j (t) 
for large t [8] . 

We remark that the exponential growth is generally 
valid only for infinitesimal perturbation. When the per- 
turbation reaches the typical size of the system, the dif- 
ference between the perturbed and the unperturbed tra- 
jectory tends to saturate. Thus, for finite amplitude per- 
turbations the mean response is the average of terms that 
remain of order O(l), and less statistics is required to 
obtain convergence. In this sense the mean response to 
finite perturbation is more representative of the behavior 
of a single perturbation than in the infinitesimal case. 

On the other hand, even if equation (13) is formally 
valid for arbitrary large perturbations, for practical use 
an upper limit exist due to finiteness of statistics. To 
predict the relaxation of a perturbation 5x(0), one needs 
sufficient statistics for the convergence of p(x(0) — Sx(0)). 
This request is more severe in systems where large fluc- 
tuations are suppressed. An example is provided by the 
stochastic model (17) with 

U(x) = \x* + ix 4 (27) 
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Here the pdf has sub-Gaussian tails, and we observe the 
opposite behavior of the system (19), as shown in Fig. 5. 
While in the case with exponential tails we have a good 
statistical convergence for a perturbation greater than 
2(7 in the second system this perturbation is too large to 
obtain convergence even with huge statistics (10 9 runs). 
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FIG. 1. Correlation function of the z variable of Lorenz 
model (solid line) compared with the mean response to differ- 
ent perturbations of the same variable. Szo = 10~ 2 ct (dashed 
line), Sz — a (dotted line), with a = \J (z 2 ) — (z) 2 — 8.67. 
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FIG. 2. Mean response of the stochastic differential equa- 
tion dx/dt = ~B(x) + s/2D£{t), with D = 1, B(x) = x for 
\x\ < 1 and B(x) = 1 for \x\ > 1, to different perturbations: 
large 8x — 2.3a (+) and infinitesimal 8x — 7.6 x W~ [i a 
(x). In both cases the mean response is exactly predicted by 
the correlator < x(t)f(x(0)) > (dashed line for Sxo = 2.3(7 
and dotted line for Sxo = 7.6 x 10~ 3 ) according to Eq.(13) 
while the simple correlation {x(t)x(0))/a 2 (solid line) just 
gives an estimate of the relaxation time. In the inset we show 
the invariant probability distribution p(x)a versus x/a with 
(7 = \J (x 2 ) — (x) 2 = 1.32. Statistics is over 10 6 independent 
runs. 
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FIG. 3. Mean response R(t) = (SE(t) /SE(0)) of the total 
energy E(t) = l/2^]|u n | 2 of the shell model (22) to dif- 
ferent amplitude perturbations: 8E(0) = 5.5 x 10~ 3 (+), 
SE(0) = 1.7 x 10" 3 (x), 5E(0) = 4.5 x 10 -4 (*). Varying 
the amplitude of the initial perturbation different relaxation 
rates are observed, and the response function is roughly sim- 
ilar to the correlation function of the corresponding largest 
perturbed shell: shell n — 12 (solid line), shell n = 14 (dashed 
line) , shell n = 16 (dotted line) . 
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FIG. 4. Halving times Ti/ 2 of the mean response 
to different amplitude perturbations of the total en- 
ergy E(t) = l/2Y / \un\ 2 of the shell model (22): 
R(t) = (SE(t)/SE(0)). Solid line represents the dimensional 
scaling Ti/ 2 ~ SE. 
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FIG. 5. Mean response of the stochastic differential equa- 
tion dx/dt = -B(x) + V2D£(t), with D = 1, B{x) =x + x :i , 
to different perturbations: finite Sxo = 1.5a (+) and infinites- 
imal 5xq = 1.5 x 10 _2 a (x). The mean response is exactly 
predicted by the correlator < x(t)f(x(0)) > (dashed line for 
Sxo = 1.5a and dotted line for Sxo = 1.5 x lCF 2 er) according 
to Eq.(13), while the simple correlation (x(t)x(0)) /a 2 (solid 
line) just gives an estimate of the relaxation time. In the inset 
we show the invariant probability distribution p(x)a versus 
x/a with a = \j (x 2 ) — (x) 2 = 0.68. Statistics is over 10 6 
independent runs. 
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